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Various definitions for the QCD Debye mass and its evaluation are reviewed in a 
non-perturabtive framework for the study of screening of general static sources. 
While it is possible to perturbatively integrate over scales ~ T and thus construct 
a 3d effective theory, the softer scales ~ gT and ~ g 2 T are strongly coupled 
for temperatures < 10 7 GeV and require lattice simulations. Within the effective 
theory, a lattice treatment of screening at finite quark densities p.<4/T is also 
possible. 



1 Introduction 

One goal of heavy ion collision experiments is to probe the deconfined phase 
of QCD, in which matter is expected to exist in a state of quark-gluon plasma. 
At present, we are still far from a satisfactory understanding of its properties, 
and hence of unambiguous signals for its detection. Qualitative expectations 
are based on asymptotic freedom, i.e. the coupling weakening with increasing 
temperature. However, for any reasonable temperatures below 10 11 GeV, per- 
turbative treatments are badly convergent due to infrared-sensitive contribu- 
tions. On the other hand, the only non-perturbative first principles approach, 
four-dimensional (4d) lattice simulations, is severely restricted in the presence 
of light fcrmions, and as of yet fails completely at finite baryon density. 

Considerable progress has been made over the last years in studying the 
screening of static sources in an equilibrated plasma by simulating dimension- 
ally reduced QCD, as well as 4d pure gauge theory. In this contribution these 
developments are reviewed with special emphasis on Debye screening, which 
is believed to play a central role for deconfinement. While at present it is not 
clear whether screening of static charges can be directly probed experimen- 
tally, it is a well-defined physical plasma property, providing the dynamical 
length scales over which interactions are effective. 

One difficulty with Debye screening in QCD is properly defining it. After 
briefly recalling the QED analogue and pointing out the differences in QCD, 
a non-perturbative framework to study screening of a general static source 
is given. Next, an effective theory description for hot QCD is reviewed, per- 
mitting a convenient non-perturbative study and interpretation of screening 
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lenghts, before returning to the question of a gauge-invariant definition and 
non-perturbative evaluation of Debye screening. Finally, Debye screening in 
the presence of a finite baryon density is discussed. 



1.1 Debye screening in QED 

Consider first a QED plasma in equilibrium, into which an external static 
charge is introduced. The modified Hamiltonian is then 

H' = H + H ext = i ((E + E ext ) 2 + B 2 ) . (1) 



In linear response theory the perturbation of the electric field is 

ft 
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S{Ek(y)) = Urn i / dt' Tr (p[R cxt (x' , t'), Ei(x, t)]) 
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lim -i fdt' f d 3 x'Ef t (x')([E i (x,t),E j (x',t')]) . (2) 
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With Ei(x) = — V(V(r), the potential of the static source is determined by 
the correlation of electric fields, or equivalently by the Fourier transform of 
the Aq propagator, 

. . n f d 3 k e ikr e - mDr 

V ^ = QJ (2^ fc2 + n oo( o,k) =Q ^r- (3) 

Thus, the electric field of the source is screened by the charges in the plasma, 
and the Debye mass, or inverse screening length, is defined by the pole of the 
photon propagator, 

n o(A:oO,k 2 = -m 2 D ) =m 2 D . (4) 

The leading order one-loop diagram gives m° D — eT/\fi + 0(e 2 T). On the 
other hand, there is no magnetic screening to all orders of perturbation theory, 
i.e. Uiiiko = 0,k^ 0) = 0. 

1.2 Non-abelian Debye screening 

Following the same steps in an SU(N) theory with Nf fermions and evaluating 
the pole of the one-loop Aq propagator, one finds the leading order result 

e _m D r „ fN NA 1/2 
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However, at next-to-leading order the problem becomes non-perturbative. 
The general form of the series in g can be shown to be □ 

™d = m° D + ^9 2 T In ^| + c N g 2 T + 0(g 3 T) , (6) 

which is non-analytic in the coupling constant. While the coefficient of the 
logarithm is fixed perturbatively, cat is entirely non-perturbative. The reason 
is that, starting from this order in g, the non-abelian Aq couples to the soft 
magnetic gluons Ai ~ g 2 T, and hence becomes sensitive to the non-abelian 
infrared divergencies in the magnetic sector for which there is no perturbative 
cureu: Ha(k = 0, k — > 0) ~ g 2 T ^ 0, with contributions from all loop orders 
two and larger. 

This raises the conceptual problem whether a perturbative definition of 
Debye screening in QCD is at all sensible. Moreover, neither IIoo nor the 
colour-electric field E are gauge-invariant physical concepts in a non-abelian 
theory. It can be shown that the pole of the Aq propagator is gauge-invariant 
order by order in a (resummed) perturbation series u. However, this does not 
guarantee the existence of a pole in the full propagator. Even if it exists, its 
non-perturbative relation to physical quantities remains unclear in a situation 
where the potential of a single static charge is not defined. 

There are then two ways to deal with this problem non-perturbatively. 
One is to assume a pole in the full propagator and to determine its value 
from the exponential fall-off of the gauge-fixed propagator in space. Resulis 
of this approach are presented in another contribution to this conference 
The second option, which will be discussed here, is to seek a manifestly gauge- 
invariant definition of Debye screening. 

1.3 Non-perturbative description of screening 

Let us first recall the general definitions needed to non-perturbatively describe 
screening of any gauge- invariant source. Consider gauge- invariant, local op- 
erators A(x). Static equilibrium physics is described by spatial correlation 
functions of euclidean time averages 

C(|x|) = (A(x)A(0)) c ~ e~ M l x l, A(x) = T drA(x, -ir) . (7) 

Jo 

These fall off exponentially with distance. The "screening masses" M have 
a precise non-perturbative definition: they are the eigenvalues of the space- 
wise tranfer matrix in the corresponding lattice field theory. Physically, they 
correspond to the inverse length scale over which the equilibrated plasma is 
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sensitive to the insertion of a static source carrying the quantum numbers of 
A. Beyond 1/M, the source is screened and the plasma appears undisturbed. 
With these definitions, all screening lengths corresponding to gauge-invariant 
sources can be computed on the lattice in principle. 

2 Effective theory description 

If we are interested in length scales larger than the inverse temperature, 
|x| ~ 1/gT ^> 1/T, the situation simplifies considerably. In this case the 
integration range for euclidean time averaging in Eq. (Q) becomes very small, 
and the problem effectively three-dimensional (3d). The calculation of the 
correlation function C(|x|) can be factorized: the time averaging may be per- 
formed perturbatively by expanding in powers of the scale ratio gT/T ~ g, 
which amounts to integrating out all modes with momenta ~ T and larger 
(i.e. the non-zero Matsubara modes, in particular the fermions). This proce- 
dure is known as dimensional reduction The correlator C(|x|) of 3d fields 
is then to be evaluated with a 3d purely bosonic effective action, describing 
the non-perturbative physics of the modes ~ gT and softer, which can be 
done on the lattice. The screening masses M now are the eigenvalues of the 
3d tranfer matrix. The 3d purely bosonic simulations are much easier and 
thus more accurate than simulating directly in 4d. However, the perturbative 
integration over hard modes and neglecting higher dimension operators in the 
effective action introduces a relative error, which at two-loop level is El 

§ ~ Otf) . (8) 

In ±he treatment of the electroweak phase transition, this error is less than 
5933, but for hot QCD the coupling and thus the error is larger. An estimate 
of its size will be given later, when results from the full and the effective theory 
are compared. 

The effective theory emerging from hot QCD by dimensional reduction is 
the SU(3) adjoint Higgs model with the action 

S = J d 3 x jiTr(FijFij) + Tr(DiA ) 2 + (m^T^A 2 ,) + A 3 (Tr(Ao) 2 | . (9) 

Apart from the SU(3) gauge symmetry, the action respects SO (2) planar rota- 
tions, two-dimensional parity P, charge conjugation C and A$ -reflections R. 
Its parameters are via dimensional reduction functions of the four-dimensional 
gauge coupling g 2 , the number of colours N and flavours Nf, the fermion 
masses and the temperature T. In all of the following fermion masses are 
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Figure 1. Spectrum of screening masses in various quantum number channels at low and 
high temperatures. Filled symbols denote 3d glueball states. 



assumed to be zero, but in principle any other values may be considered as 
well. At leading order in the reduction step one then has 

gl=g 2 (fi)T, X 3 = -L s (6 + N-N f )g*(ji)T. (10) 

The reduction step has been performed to two-loop order i at which param- 
eters have relative accuracy 0{g A ). Specifying T and a renormalization scale 
p. completely fixes the effective theory. 

If we are interested in the longest scales of order ~ l/g 2 T, one may be 
inclined to integrate out the fields Aq ~ gT as well, leaving a 3d pure gauge 
theory. In order to decide which is the correct effective theory in practice, 
simulation results for screening lengths will now be compared between the 
reduced and the full theory. _ 

Results from simulations E3 of the gauge-invariant screening spectrum 
with the effective theory Eq. (^) are displayed in Fig. ^. Screening masses are 
classified according to the symmetires of the 3d Hamiltonian by the quantum 
numbers Jjj . The open symbols show states dominated by scalar operators 
such as Tr(Ap), Tr(AoF^ 2 ) etc., whereas full symbols show states receiving 
only gluonic contributions Tr(F^ 2 ) Mp-ift * s a remarkable finding of detailed 
mixing analyses in several models EiHl3£°l, that the latter are quantitatively 
consistent with the 3d glueball states t3 and completely insensitive to the 
presence of the Aq scalar field. Because of this pronounced non-mixing, we 
can state a first important result: For any reasonable temperatures the largest 
correlation length of gauge- invariant, local operators belongs to the Aq degrees 
of freedom and not to the A+, in contrast to the naive parametric picture. This 
demonstrates that the physics from the scale gT down to g 2 T is completely 
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Figure 2. Comparison of screening masses inrhpt SU(2) (left) and SUte)i|(*ight) pure gauge 
theory as determined in 4d (empty symbols) £3 an d 3d (full symbols) oo effective theory 
simulations. 



non-perturbative, and hence Aq may not be integrated out perturbatively. 
Only for temperatures T > 10 7 T C becomes the coupling small enough that the 
parametrically expected ordering of screening masses is realized. 

Next, we are interested in the accuracy of the reduced theory. Fig. || 
compares the results for hot SU(2) gauge theory as obtained in the full lattice 
theory Ej with those from the effective theory Ej. Note that the effective theory 
is only valid up to its cut-off M/T ~ 2tt and above this level disagreement 
is to be expected. The apparent excellent agreement for the two lowest lying 
states should be taken with care given some systematic uncertainties. In 4d 
simulations, the projection properties of operators and the scaling behaviour 
are more difficult to control than in 3d. In the effective theory, there is 
an ambiguity in the choice of renormalization scale. The uncertainty in the 
comparison may thus be as large as 20%. For the case of SU(3), one finds 
again quantitative agreement in the largest correlation length but about 20% 
deviation in the next shorter one. 

One advantage of the effective theory approach is its easy inclusion of 
fermions, which are treated analytically and just amount to changing the pa- 
rameters Eq. ( |l0| ) of the action to be simulated. However, the screening masses 
M/T are found to increase with Nf EL3, and at Nj — 4 the lowest of them 
is already of the order ~ X/lixT and hence not expected to be quantitatively 
accurate. Moreover, close to T c the fermionic modes in the 4d theory feel 
the chiral phase transition, such that for T < 2 — 3T C the correlation length 
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obtained from pion operators becomes the longest onet-2. In this situation 
fermions may not be integrated out and the purely bosonic effective theory is 
invalid. We may thus conclude that dimensional reduction gives a reasonable 
description of the largest correlation lengths of hot QCD for temperatures 
> 2 — 3T C , with an error < 20% that decreases with temperature. 



3 The Polyakov line and the heavy quark potential 

A gauge- invariant, non-local operator that has been suggested in the literature 
to describe Debye screening is the Polyakov line 

1 

L(x) = -TrPexp-ig J clrA„(x,r) . (11) 

Its spatial correlator is related to the free energy of a quark anti-quark pair 
in the plasma, and thus to the static potential at finite temperature E£l 

(LHr)L(O)) V(r) 

= exp -PF q g = exp — . (12) 

The idea is that by considering a source together with an anti-source gives 
a manifestly gauge-invariant operator related to heavy quark systems. In 
leading order perturbation theory one has L = 1 — g 2 /2Tr(J drAo) 2 + . . ., 
and hence the correlator at this order is given by the exchange of two electric 
gluons leading to 

(N 2 - Do 4 c~ 2m D r 

with the leading order Debye mass m° D . However, this state of affairs is spoilt 
by higher loop diagrams, as several authors have pointed out EZrE3. Fig. ^| left 
shows a diagram of order g 8 , where now Aq can couple to two magnetic gauge 
fields Ai , and hence the exponential fall-off of this diagram is governed purely 
by magnetic gauge fields. Indeed, through a fermion loop a similar coupling is 
possible even in QED, cf. Fig. fright. There, magnetic gauge fields are strictly 
massless, and the correlator falls off algebraically rather than exponentially. 
It thus fails to produce the Debye mass for electric fields even in QED. 

With our previous discussion of screening masses at hand, it is not difficult 
to understand how the Polyakov line correlator should fall off away from T c . 
L contains a sum over all J PC sectors in the fundamental representation. 
Furthermore, because of gauge-invariance it couples only to gauge-invariant 
intermediate states. Its correlator should therefore asymptotically fall off with 
the lightest gauge-invariant screening mass of the spectrum. Indeed, a recent 
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Figure 3. Higher loop diagrams contributing to the Polyakov line correlator in QCD (left) 
and QED (right). 



4d simulation of the SU(3) Polyakov line correlator E2I finds it exponential 
decay to be governed by M/T s» 2.5(3), which is fully compatible with the 
lightest 0+ + mass given in the table in Fig. 0. The same observation is made 
to rather high precision in a hot (2+l)-dimcnsional gauge theory Ell. 



4 A gauge- invariant definition for mjj 



A few years ago it was suggested to use a symmetry to-distinguish between 
the electric and the magnetic sectors of the theory Ea, namely euclidean 
time reflection T : t — > — t. This operation changes the sign of the elec- 
tric gauge fields, Aq — > — Aq, but not that of the magnetic ones. The 
non-perturabtive Debye screening length 1/rno may then be defined to be 
the largest correlation length measured from a correlator (A(x)B(O)), where 
A, B are gauge-invariant operators which are local in 3-space and odd un- 
der T. Examples of four-dimensional operators with these properties are 
Tr(ImL), Tr(ImLF 12 ), Tr^F^), Tr(F 03 F 12 ) etc. 

In the 3d effective theory, where "time" has been reduced away, the rem- 
nant of euclidean time reflection is the scalar reflection symmetry R : Aq — » 
— Aq. The Debye mass then corresponds to the lightest eigenstate with quan- 
tum number R = — . It is now easy to go back to the results for the spectrum 
in Fig. |l| and identify the lowest mass in the 0Z + channel as the Debye mass 
according to this definition^. The corresponding operator is Tr(ArjFij), which 
may be interpreted as an electric gauge field dressed by magnetic gluons to 
make it gauge-invariant. Despite its coupling to A%, i?-symmetry prohibits 
its correlation to be governed by the 3d pure gauge sector. This definition 
certainly refers to a physical quantity. However, as it couples to colour sin- 
glet sources, its relation to screening of colour flux and to deconfinemcnt of a 
quarkonium system is not obvious. 

a Alternatively, it was proposed to measure matrix elements betjween _R-odd operators and 
a Zjv domain wall, introduced by a twist in the effective action 
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Table 1. The different contributions to the Debye mass, Eq. (|3|). (1+2) refers to the first 
two terms in units of g 2 T. 





m D /g 2 T 


(1 + 2) 


C3 


0(g*T)/{g-T) 


T = 2T C 
T = 10 n T c 


1.82 (3) 

3.83 (9) 


0.514 
2.165 


1.65(6) 
1.65(6) 


-0.35(8) 
0.015(90) 



4-1 Perturbative and non-perturbative contributions to 

We may now ask to what extent this Debye mass can be computed in pertur- 
bation theory. First, note that the perturbative expansion follows the same 
pattern as in Eq. (Q). At tree level, the correlator is given by exchange of 
one electric and two magnetic gluons. Since the latter are massless at tree 
level, the leading order fall-off is yet again m° D . However, at higher orders the 
contributions are different from those to the A$ propagator. _ . 

There is a recipe to compute the coefficient cat non-perturbatively 113: if 
temperature is sent to infinity, the Aq becomes infinitely heavy and decouples, 
so it may be integrated out. In this case, Aq turns into a static external source, 
and the only dynamical fields left are the Ai of a 3d pure gauge theory. We 
thus have 

-m^lxl 

(Tr(A o F ij )(x)Tr(A o F lm )(0)) — > _^(Fg(x)W ab (x, O)F^(O)) , (14) 

where W a b{~x-, 0) is an adjoint representation Wilson line. 

As a result, our correlator factorizes into the free Aq propagator and a 
correlator of the magnetic field strength, connected by an adjoint Wilson line 
representing the static Aq . The field strength correlator falls off exponentially 
~ exp — (A + c/v)|x|, where A contains a logarithmic -divergence which gets 
absorbed by the renormalisation of the Aq self-energy 113. 

Measuring cn in a 3d pure gauge theory 113 and comparing with the mea- 
surements of the full mj) we can now summarise the various contributions to 
the SU(3) Debye mass in Table |j We observe that at reasonable temperatures 
the next-to-leading order coefficient C3 is larger than the leading contribution, 
and thus the Debye mass is entirely non-perturbative. Again the reason is its 
coupling to the non-perturbative magnetic sector ~ g 2 T. Note also that the 
0(g 3 T) contributions are small compared to the g 2 T contributions. Together 
with the previous observations about scales involved, one may then conclude 
that the physics of screening for all temperatures of pratical interest is domi- 
nated by the scale ~ g 2 T. Only at asymtotically high temperatures T > 10 7 T C 
is the perturbative picture restored. 
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5 Debye screening at finite baryon density 



In heavy ion collisions the initial state has non-vanishing baryon number, 
and so does any subsequent plasma state. We would therefore like to know 
if and how a chemical potential for quarks affects the physics of screening. 
Unfortunately, finite baryon density cannot be addressed by standard lattice 
QCD because of the so-called "sign-problem" : when a chemical potential term 
is added to the theory, the fermion determinant becomes complex prohibiting 
Monte Carlo importance sampling, which requires a strictly-positive measure. 
To date, no working cure has been found for this problem IH3. 

Recently it was demonstrated that in dimcnsionally reduced finite density 
QCD the sign problem is in fact numerically tractable c3. Inclusion of a 
chemical potential term for quarks leads to one extra term in the action Eq. (Q) 
and changes the tree level Debye mass, 



S ->S- 



.H Nf 



d^xTrAS 



1 



fi \ 2 3N f 
ttT/ 2N + N f 



(15) 



This action is complex and cannot be simulated with standard methods. In- 
stead, a reweighting procedure has to be performed, where the complex term 
of the action is absorbed into the observable, so that 



(O) 



(Oe~ 



cos(S" M ))o 



(16) 



and the averaging is done with fi = 0. The sign problem, i.e. cancelling 
contributions to the expectation value, occur whenever 1. Fortunately, 

the width of the distribution of this quantity, growing as ~ (/i/T)^ 1 / 2 , is 
narrow enough to permit large enough volumes-for simulated masses to reach 
their infinite volume levels, as long as [i < AT 113. 

There is one important change of the symmetries of the theory, which 
is relevant for Debye screening L3: the chemical potential term violates R 
and C, and thus only P is left as a good symmetry to classify the screening 
states. Since the definition of mo was based on the use of R, it seems no 
longer applicable in the finite density situation. In the 3d effective theory, 
the corresponding state can mix with any lighter R — + states and decay. 
However, for fj, = 0, we found the Debye mass to be the lowest state in the 
0I + channel. This is also the lowest state with P = — . As Fig. ^ shows, this 
remains to be the case as fi/T is switched on, and by continuity one may now 
identify the lowest 0~ mass to correspond to the Debye mass. 

Finally, it has been demonstrated, that screening masses for fi < T should 
be accessible to 4d simulations with imaginary fi and analytic continuations. 
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Figure 4. 

switched on. The triangles in the right panel give the Debye mass according to the gauge 
invariant definition in the last section. 



6 Conclusions 

The analogy of Debye screening in QED and QCD is rather limited, because in 
QCD it is not possible to define the electric field induced by an isolated charge: 
neither of these are gauge-invariant physical concepts. The problem then is 
to specify what we mean by Debye screening in QCD. Depending on the def- 
inition, one gets different answers and different interpretations of the results. 
Defining the Debye mass through the pole of the full A propagator has the 
disadvantage that its existence is not guaranteed, and its non-perturbative 
connection to physical quantities is not clear. Employing the Polyakov line 
correlator describing a heavy quark anti-quark pair does not probe Debye 
screening: it mixes electric and magnetic sectors, failing to fall off with nip 
even in QED. A gauge- invariant definition is possible through correlations of 
operators odd under euclidean time reflection, which reflects Ao but not Ai. 

In general, the screening of static sources in equilibrium is a well-defined 
problem that has been studied non-perturbatively by lattice QCD and with 
effective field theory methods. A detailed picture of the static length scales for 
non-abelian plasma physics has emerged. An important lesson for construct- 
ing effective theories is that, for all temperatures T < 10 7 GeV, the scales 
<~ gT and ~ g 2 T are not separable non-perturbatively. All static screening 
physics seems to be dominated by the dynamics on soft scales ~ g 2 T. 

Finally, dimensionally reduced QCD gives a realistic description of 
static phenomena down to temperatures of a few T c , and is able to non- 
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perturbatively accomodate chemical potentials for quarks \i < AT. It can 
therefore address the physics of screening in the parameter range relevant 
for heavy ion collisions. 
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